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Abstract- 

We  have  devised  a  system  to  automatically  detect,  manage 
and  track  contacts  with  a  stationary  line  array.  The  array  has  a 
large  number  of  elements,  and  the  ability  to  process  any  subset  of 
these  elements  (up  to  a  fixed  number)  at  any  given  time. 
Detection  is  accomplished  by  fixed  detection  apertures  spaced 
along  the  array  such  that  the  apertures  create  a  “tripwire”  range 
parallel  to  the  length  of  the  array.  Contacts  crossing  the  tripwire 
range  are  handed  off  to  an  automated  tracking  aperture  which 
adjusts  its  aperture  size  and  aperture  center  along  the  array  to 
track  the  contact  using  the  best  subset  of  elements.  The  adaptive 
detection  scheme  conserves  limited  tracker  resources  by  limiting 
the  trackers  to  those  contacts  which  appear  to  be  heading 
towards  the  line  array.  This  system  has  beenO  implemented  on  a 
Mercury  Adapdev  system  and  has  been  tested  using  simulation  of 
multiple  approaching  contacts  with  interfering  distant  contacts. 

I.  Introduction 

We  have  devised  a  system  to  automatically  detect,  manage 
and  track  contacts  with  a  stationary  line  array.  The  array  has  a 
large  number  of  elements,  and  the  ability  to  process  any  subset 
of  these  elements  (up  to  a  fixed  number)  at  any  given  time. 

The  detection  scheme  is  designed  to  conserve  limited  tracker 
resources  by  limiting  the  trackers  to  contacts  that  appear  to  be 
heading  towards  the  line  array.  Detection  is  accomplished  by 
automated  broadband  energy  followers  operating  on  fixed 
detection  apertures  spaced  along  the  array  such  that  the 
apertures  create  a  “tripwire”  range  parallel  to  the  length  of  the 
array.  A  detection  management  system  was  built  to  handle  the 
detection  aperture  contacts.  Contact  followers  (CF)  are 
automatically  assigned  to  all  detections  in  these  apertures.  The 
detection  system  can  track  numerous  contacts  in  bearing  and 
range  in  both  near  field  and  far  field  in  each  detection  aperture. 

Contacts  that  cross  into  the  tripwire  range  and  are  closing  are 
handed  off  to  a  tracking  aperture,  which  adjusts  its  location 
along  the  array  to  track  the  contact  using  the  best  subset  of 
elements.  Candidate  CF’s  are  used  as  detections  to 
automatically  assign  the  range-focused  trackers.  The  tracker 
aperture  is  then  allowed  to  “rove”  along  the  length  of  the  array 
to  keep  up  with  contact  dynamics.  The  tracker  also  adjusts  its 
aperture  size  with  the  changing  range  of  the  contact.  The 
tracker  incorporates  its  own  range-focused  beamformer 
utilizing  a  three  beam  interpolation  in  azimuth  and  inverse- 
range  space.  It  uses  an  adaptive  gain  alpha-beta  filter  for  state 
estimates.  The  predicted  state  estimates  provide  range  and 
bearing  estimates  to  reset  the  tracker  beamformer  aperture  to 
the  new  position.  There  are  a  limited  number  of  these  trackers 
since  they  require  front-end  beamformer  functionality  and  run 
at  a  faster  rate  to  handle  contact  dynamics.  The  track 


management  function  oversees  the  tracker,  and  adds  and  drops 
them  based  on  multiple  criteria  such  as  contact  opening  or 
closing  and  signal  to  noise  information. 

This  system  is  implemented  on  a  Mercury  Adapdev  system 
and  has  been  tested  using  simulation  of  multiple  approaching 
contacts  with  interfering  distant  contacts. 

This  paper  focuses  on  the  tripwire  detection  and  localization 
scheme.  In  Section  II,  the  range  focused  tripwire  concept  is 
developed  in  detail.  Section  III  discusses  the  application  of  the 
tripwire  to  the  long  line  array  detection  system.  Section  IV 
shows  simulation  results. 

II.  Range  Focused  Tripwire  Processing 

A.  Range  Focused  Beamforming 

A  key  concept  in  the  definition  of  this  detection  scheme  is 
design  of  the  focus  range  sets  for  a  long  line  array,  specifically 
line  arrays  of  sufficient  length  such  that  the  near  field  extends 
out  past  ranges  of  practical  importance 

iLl  >r 

^  IK  significan  t 

Where  L  is  the  aperture  length;  f  is  frequency  and  C  is  the 
speed  of  sound  in  water  as  defined  in  [1]. 

In  order  to  provide  full  coverage  of  ranges  of  practical 
importance  that  would  be  degraded  by  a  plane  wave  (infinite 
range)  beam  set,  one  or  more  beam  sets  focused  at  closer  finite 
ranges  will  be  required.  Define  a  focus  range  set  R  /  (6)  as 

specifying  a  set  of  ranges  as  a  function  of  bearing  0  referenced 
to  the  center  of  a  line  array.  Let  j  =  {0,  1,2,  ...,  N}  be  integers 
representing  nested  consecutive  focus  range  sets  such  that 

Rq{9)> Rx(e)> R2{e)>....> r  .(e)  \/e-.  o°  <0<iso° 

The  focus  range  set  If  (9)  ,  j=  {0,  1,2,  . . .  N}  represent  range 
and  bearing  coverage  with  respect  to  a  line  array.  The  farthest 
range  as  a  function  of  bearing  coverage  is  specified  by  R{]  (  9) 
and  the  closest  range  as  a  function  of  bearing  coverage  is 
specified  by  RN  (i 9 ) . 

We  begin  our  study  of  focus  range  design  by  examining  the 
time  delay  equation  between  a  point  in  space  (bearing  0  and 
range  R)  and  an  array  element.  This  analysis  will  yield  insight 
into  focus  range  design.  [2],  [3],  and  [4]  contain  additional 
treatments  of  range  focusing. 

We  will  assume  a  long  line  array  centered  on  the  x-axis.  The 
point  in  space  which  we  wish  to  beamform  to  is  assumed  to  be 
at  a  range  R  from  the  center  of  the  line  array  and  at  an  angle  0 
from  the  center  of  the  line  array.  Let  d;  be  the  ith  sensor 
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location  on  the  line  array.  Fig.  1  illustrates  this  geometry. 
Using  the  law  of  cosines  we  can  write  the  following  equation 
representing  the  coherent  contact  signal’s  time  of  arrival 
difference  x,  between  sensor  d;  and  the  center  of  the  line  array, 
given  speed  of  sound  C. 

,  .  R-Jr2  -  2-R-d:  -  cos O+df 

Ti(0’R)  = - 1 - 7TJ - -  CD 


The  first  term  is  the  standard  delay  for  plane  wave 
beamforming  of  a  far  field  contact  at  angle  0  to  the  center  of 
the  line  array.  The  second  term  is  the  time  delay  correction  for 
a  near  field  contact  at  range  R  from  the  line  array.  The  second 
tenn  is  of  primary  interest.  This  term  is  effectively  added  to  a 
plane  wave  beamformer  to  produce  a  focus  beam  at  range  R 
and  bearing  0. 


Figure  1:  Near  field  geometry 

The  following  trick  results  in  a  very  accurate  approximation 
and  yields  insight  into  focused  beamforming.  For  the  purposes 
of  this  study,  and  in  order  to  support  the  two  approximations 
to  follow,  we  shall  assume  that  R  »  d  ,•  •  cos  #  ,  that  is,  the 
range  is  much  larger  than  the  array  element  spacing. 


A  little  algebra,  completing  the  square,  and  using  of  the 

9  9 

identity  sin”  #  +  cos  9  =  1  yields  the  following 
R  -  tJ(R-  dj  ■  cos#)2  +  df  ■  sin2  9 


Ti(9,R)  =  - 


C 


(2) 


Factoring  out  the  dominant  tenn  yields 
R  —  (R  —  dj-cos9)-_  1  + 


-(M=- 


1 


df  ■  sin2  9 


(R  -  dj  ■  cos  9) 


C 


Now  we  apply  the  following  approximation  to  (3) 

Vl  +  a  ~\+—  for  a  « 1 
2 

Which  yields 


(3) 

(4) 


R  -  (R  -  dj  •  cos 0)  ■ 


ZiM- 


df  ■  sin2  0  ''' 

2  \R-dj-  cos#)'  j 


C 


(5) 


and  finally  after  some  algebra 

Tl(e,*)~drmse 


df  ■  sin”  0 


C  2  (R-dj  ■  cos  #)  •  C 
Equation  (6)  can  be  approximated  as 


dj-  cos  9  df  sin2  0 


C 


2-R-C 


(6) 


(7) 


We  will  now  look  at  the  time  delay  difference  between  a  beam 
in  focus  range  j  focused  to  range  R  j  (9)  and  bearing  0  and  a 

focus  range  j+1  focused  to  range  Rj+l  (9)  and  bearing  0. 


i ^rJ-MRj)-- 


df 

2-C 


-•si: 


itr  6- 


Rj{e)  Rj+X{e) 


(8) 


Now  if  we  require  the  scalloping  loss,  i.e.,  the  maximum 
reduction  in  power  between  any  pair  of  beams  focused  on  the 
same  bearing  in  adjacent  focus  range  sets,  to  be  constant  for 
all  bearings  9  and  all  focus  range  sets  j  ,  then  we  require 


sin2#- 


Rj{e)  RjJe) 


=  Const,  V#Vjr  (9) 


Let  us  define  a  constant  parameter  A  as  shown  below: 

1  1 

A  =  - r - v - r - v  (10) 

R/+1(90°)  R  j  (90  0  ) 

The  parameter  A  specifies  the  desired  spacing  between  focus 
range  sets  at  90°  to  be  a  constant  A  and  thus  effectively 
specifies  the  desired  range  scalloping  loss  between  adjacent 
focus  ranges  sets. 

We  can  now  define  the  relationship  between  the  focus  ranges 
of  each  pair  of  beams  in  adjacent  focus  range  sets  as  shown 
below,  resulting  in  a  fixed  range  scalloping  loss  with  respect 
to  focus  range  set  R  /+1  (#)  for  all  common  bearing  angles  0. 


1  _  1 

Rj+xiefTjWl 


- - - for  0°  <#<180° 

sin2(#) 


(ID 


B.  One  Degree  of  Freedom 

Given  a  particular  focus  range  set  R ;  (#) ,  the  following 
relationship  with  the  other  focus  range  sets  will  guarantee  a 
fixed  range  scalloping  loss  at  common  bearing  angle  0 
between  all  adjacent  focus  range  sets  at  all  bearing  angles. 
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1 


1 


1 


A 


"RjM  RjiP)  Rfo)  Rj-M"  sin2(f?)  (12) 
for  0°  <6<  180° 


For  convenience  we  will  define  the  range  R  !n:  where  we  want 

to  have  a  constant  power  response  in  tenns  of  the  parameter  A 
and  a  selectable  parameter  b. 


There  is  no  implied  restriction  on  the  definition  of  focus  range 
set  R  •  (O) ,  except Rj(d) can  not  equal  0.  We  will  allow  focus 

ranges  to  become  large  negative  ranges  (beyond  infinite  range) 
for  now.  Looking  at  (12)  it  becomes  apparent  that  we  can 
arbitrarily  define  the  range  versus  bearing  of  one  focus  range 
set  and  then  under  the  assumption  of  equal  scalloping  loss 
everywhere  the  other  focus  range  sets  will  be  detennined. 

Thus  we  claim  that  the  definitions  of  focus  range  sets  with  the 
property  of  maintaining  the  same  fixed  scalloping  loss 
between  adjacent  focus  range  sets  at  common  bearing  angles 
has  one  degree  of  freedom,  i.e.,  the  definition  of  one  of  the 
focus  range  sets.  In  theory  there  are  an  infinite  number  of 
possible  definitions  of  focus  range  sets.  We  shall  concern 
ourselves  with  the  important  class  of  focus  range  set  definition 
where  we  maintain  a  constant  power  response  at  a  given  range 

R  jh  for  all  bearing  angles  0.  (Note  that  R  fix  does  not  have  to 

be  a  member  of  any  of  the  defined  focus  range  sets  R  • .)  The 

natural  algorithm  of  forming  a  plane  wave  beamformer  at 
infinite  range  followed  by  a  number  of  closer  finite  focus 
ranges  is  a  subset  of  this  class. 


~~  =b  -  A  (15) 

Rfix 


We  will  now  define  where  we  want  to  place  R  fix  with  respect 
to  our  focus  range  sets  j=0,  1,2,  . . .,  N  We  will  place  R  /ix  at  a 
range  set  index  of  “a”  Therefore  Ra  (9)  =  R  ;7v  and  using  (14) 
we  can  write  the  following  equation. 


R 


fix 


rM 


=  a 


If  we  now  solve  (16)  for 
1  1 


1 


*o(*) 

A 


R0  (&)  R  fix 


sin 


(16) 

we  have  the  following 
(17) 


C.  Constant  Power  Response  at  a  Fixed  Range  Across  Bearing 

We  shall  now  derive  expressions  that  define  focus  range  sets 
that  contain  a  constant  power  response  at  a  fixed  range  for  all 
bearings. 


With  (10)  in  mind  the  following  equation  satisfies  (11)  for  all 
focus  range  sets  represented  by  the  integers  j=0,  1,2,  . .  .N 
where  R0  (O)  is  the  farthest  focus  range  set. 


1  1 
Rj+i{0)~W) 

A 


(13) 


sin  2  ( 9 ) 


Let  i=0  and  thus  R0  (0)  is  the  farthest  focus  range  set. 
Equation  (13)  can  be  rewritten  as  follows. 

1  1 


R,(o)  rM 


=J 


(14) 


sin 


:(*) 


Note  that  j  is  effectively  a  linear  function  of  focus  range  set 
index.  Integer  values  of  j  correspond  to  focus  range  sets  we 
have  designed  to  have  a  fixed  scalloping  loss  between 
adjacent  sets  for  all  bearings,  j  has  the  range  of  0  >  j  >  N 
where  j  =  0  is  the  farthest  focus  range  set  and  j  =  N  is  the 
nearest  focus  range  set. 


Substituting  (17)  and  (15)  into  (14)  we  have 


— rx-b-  A  +  a - - — 

Rj{0)  sin» 

A 

sin  2  (<?) 


(18) 


Solving  for  R  .  (9)  we  have 
sin  2  (#)  1 


rM= 


j-a  +  b-  sin“(f?) 


j  =  0,1,2...,  N  (19) 


II  Observations 

For  j  >  a,  the  focus  ranges  Rj  (9)  will  decrease  as  0  moves 
away  from  90°.  For  focus  range  indices  j  <  a  the  focus  ranges 
Rj  (9)  will  increase  as  0  moves  away  from  90°.  In  fact  they 

will  achieve  infinite  range  (plane  wave)  when  sin  2  (d)  = - — . 

b 

As  0  continues  to  move  further  away  from  90°  the  ranges  will 
become  finite  and  negative.  There  is  nothing  wrong  with  this 
mathematically  or  in  practice.  The  negative  ranges  will  be  in 
focus  at  no  realizable  positive  range,  but  will  maintain  the 
correct  relationship  with  the  other  focus  range  sets.  In  practice 
a  designer  may  wish  to  limit  negative  focus  ranges  to  infinite 


Approved  for  public  release;  distribution  is  unlimited. 


range  and  abandon  the  design  relationship  with  the  adjacent 
closer  focus  range  set. 

E.  Practical  Designs 


3)  Defining  focus  ranges  for  each  focus  range  set  by  the 
following  equation  guarantees  a  constant  scalloping 
loss  at  all  bearings  between  adjacent  focus  range  sets 
but  also  specifies  a  given  focus  range  where  the 
power  response  is  constant  for  all  bearings. 


For  many  practical  designs  the  choice  of  R  fix  would  be  one 

of  the  natural  focus  range  sets  (we  would  want  one  of  our 
focus  range  sets  to  at  a  constant  range  as  a  function  of  bearing.) 
In  this  case  the  parameter  “a”  would  be  one  of  the  integers 
between  0  and  N.  A  more  interesting,  and  in  our  case,  more 
useful  choice  is  to  make  the  scalloping  range  between  two 
focus  range  sets  be  at  a  constant  range  as  a  function  of  bearing. 
This  implies  that  contacts  more  focused  in  more  distant  focus 

range  sets  would  be  at  a  range  greater  then  R  fix  independent 
of  bearing  and  contacts  more  focused  in  closer  focus  range 
sets  would  be  at  a  range  less  then  R  fix  independent  of  bearing 

and  contacts.  For  this  case  the  parameter  a  could  take  on  the 
values  {-1/2,  1/2,  1+1/2,  2+1/2,  ...,N+l/2}. 


For  many  reasons  it  may  not  be  desirable  to  set  the  farthest 
focus  range  to  infinite  range  at  90°  (Ro(90°)).  A  reasonable 
alternative  is  to  set  the  range  corresponding  to  a  bearing  of  90° 
to  have  a  power  loss  at  infinite  range  equal  to  the  scalloping 
loss.  Set  j— 0,  0=90°  and  Ro(90°)  =2/A  in  (19)  to  obtain  the 
correct  relationship  between  the  parameters  a  and  b  for  this 
case. 

b-a  =  -  (20) 

2 

Under  this  condition  the  relationship  between  the  parameter  a 
and  the  range  R  /(r  of  constant  power  can  be  obtained  by 
substituting  (20)  into  (19)  and  letting  j  =  a . 


R 


fix 


1 

A  ■  (a  +  .5) 


(21) 


Finally  the  standard  definition  of  focus  range  sets  with  the 
farthest  focus  range  set  being  set  at  infinite  range  for  all 
bearings  (plane  wave  beamformer)  can  be  obtained  from  (19) 
by  setting  a  =  0  and  b  =  0  resulting  in 

=  J  =0,l,2-„  N  (22) 

A -J 

Note  j=0  corresponds  to  infinite  range  i.e.  a  plane  wave 
beamformer. 


F.  Summary 


1) 

2) 


Selecting  the  parameter  A  specifies  the  desired 
scalloping  loss  between  adjacent  focus  ranges. 
Maintaining  the  following  relationship  between 
adjacent  focus  ranges  guarantees  a  constant 
scalloping  loss  at  all  bearings 


1 _ 1 

Rj+x{9)~~R~(e) 


- - -  for  0°  <0<18O° 

sin2(0) 


Rj{o)  = 


sin 


\e) 


j  -  a  +  b  ■  sin  2  {d ) 


j  =  0,1,2...,  N 


4)  Careful  selection  of  the  parameters  “a”  and  “b”  allow 
customizing  the  focus  range  sets  to  meet  particular  design 
goals. 


III.  Application  to  Long  Line  Array  Tripwire 
Detection 

The  range  focused  tripwire  concept  was  applied  to  the 
problem  of  automated  detection  on  a  long  line  array  with  a 
large  number  of  elements.  In  this  case,  the  detection  system  is 
required  to  sort  out  contacts  and  pass  on  to  the  track 
management  system  only  those  meeting  the  profile  of  being  in 
close  proximity  and  closing  on  the  array.  As  a  result,  two 
range  focus  beam  sets  were  chosen,  near-field  and  far-field, 
based  on  the  choice  of  a  tripwire  range  set  defining  the 
proximity  at  which  the  contact  should  be  assigned  to  a  tracker. 

Given  the  desired  tripwire  range,  Rth,  we  select  a  =1/2  to 

set  the  constant  range  tripwire  between  range  0  and  range  1. 
As  discussed  above,  these  conditions  determine  the  focus 
ranges  for  R0(0)  and  If  (O) ,  sketched  in  Fig.  2. 


\ 


> _ 

'Array 


Detection  Aperture 


Figure  2:  Tripwire  detection  scheme 

hi  our  application  the  ranges  appear  as  in  Fig.  3.  Note  that 
the  decision  to  use  fixed  far  field  ranges  instead  of  the 
negative  ranges  discussed  in  Section  II  result  in  curved 
regions  of  the  tripwire  range  near  endfire.  Practically  speaking, 
ranging  performance  is  best  for  those  beams  focused  on  angles 
within  approximately  30  degrees  of  broadside.  Thus,  a  single 
detection  aperture  provides  excellent  coverage  of  the  tripwire 
range  within  a  60  degree  wide  azimuthal  swath  as  shown  in 
Fig.  4. 
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IV.  Simulation  results 


Three  focused  ranges 


Figure  3:  Tripwire  range  sets  for  long  line  array 

The  length  of  the  array  allows  for  multiple  detection 
apertures  to  be  processed.  The  detection  apertures  are  spaced 
along  the  array  as  shown  in  Fig.  5  and  the  detection  manager 
limits  the  extent  of  the  beams  in  each  detection  aperture  fan  to 
support  identifying  the  closest  aperture  as  a  contact  passes 
inside  the  tripwire  range,  so  that  there  is  sufficient  scalloping 
loss  between  the  outer  beams  of  adjacent  apertures  at  the 
tripwire  range. 

~  60  deg  coverage 
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Figure  4:  Single  Detection  Aperture 


Two  time  cuts  of  beamformed  data  for  an  example  scenario 
in  a  single  aperture  are  shown  in  Fig.  6.  (Note  that  the 
additional,  inner,  range  window  shown  was  added  for  use  in 
the  tracker  manager,  not  discussed  in  this  paper.)  Four 
contacts  were  simulated;  two  were  stationary  (labeled  1  and  4 
in  Fig.  6)  and  outside  the  tripwire  range,  while  the  two 
dynamic  contacts  started  at  070  and  090  degrees  respectively, 
relative  to  the  aperture  center  and  crossed  over  the  line  array. 
These  contacts  are  labeled  2  and  3. 

The  first  time  cut  labeled  “Approaching  CPA”  shows  an 
early  cut  of  the  data.  All  four  contacts  are  initially  and 
correctly  focused  in  the  R0  (far)  window,  representing  the 
furthest  focus  range.  Looking  at  the  R1  (mid)  window,  notice 
that  the  two  outside  contacts,  which  are  furthest  from  the  array, 
are  more  defocused  than  the  closer,  dynamic,  contacts  which 
become  more  focused  over  time  as  they  move  closer  to  the 
array.  The  defocusing  is  most  pronounced  in  the  innermost  R2 
(near)  window. 

The  second  time  cut  labeled  “CPA”  shows  data  from  the 
same  contacts  around  the  time  when  the  two  dynamic  contacts 
pass  CPA  (closest  point  of  approach)  directly  over  the  array 
simultaneously.  Both  dynamic  contacts  are  focused  in  the  R1 
(mid)  window  for  most  of  this  time  period,  while  the  static 
contacts  remain  focused  in  the  R0  (far)  window.  The  first 
contact  passes  over  at  the  left  edge  of  the  aperture,  moving 
rapidly  through  the  beams  as  it  approaches  and  departs.  The 
second  contact  passes  directly  over  the  center  of  the  detection 
aperture,  focusing  in  the  R2  (near)  window  at  the  tune  closest 
to  CPA  and  defocusing  everywhere  as  it  reaches  CPA 
(partially  a  result  of  the  fixed  detection  aperture.) 

V.  Summary 

We  described  the  detection  processing  scheme  used  for  a 
stationary  long  line  array  designed  to  automatically  detect, 
manage  and  track  contacts.  A  set  of  range  focused  beamsets 
were  designed  to  provide  a  constant  range  tripwire  for 
automated  identification  of  contacts  approaching  the  line 
array  and  were  replicated  over  multiple  apertures  along  the 
length  of  the  array.  A  simulation  example  was  presented  to 
illustrate  the  localization  concepts  of  this  approach. 
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Figure  5:  Multiple  Apertures  Create  Tripwire  along  the 
Length  of  the  Array 
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Figure  6:  Simulation  Example  showing  multiple  range  windows  for  a  single  detection  aperture.  Contacts  1  and  4  are 
stationary  far  field  contacts.  Contacts  2  and  3  are  approaching  the  array.  Contact  2  crosses  the  array  at  CPA  to  the  left  of 
the  detection  aperture  in  the  R1  window.  Contact  3  crosses  the  array  at  CPA  in  the  center  of  the  detection  aperture  in  the 

R2  window. 
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